Scale invariant pattern in dynamically extending lattice 
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We study the implications of coupling Langmuir kinetics (LK) process with the dynamics of 
an extending lattice. The model that we consider couples dynamically extending exclusion process 
(DEEP) with the process of random attachment and detachment of particles in the bulk of the lattice. 
We explore a dynamical regime where the boundary processes of lattice extension and particle input 
compete with the bulk process of particle attachment and detachment. This competition leads 
to scale invariant density profile of particles when expressed in terms of the relative position in 
the growing lattice. It also leads to phase coexistence and shocks in the bulk of the lattice. We 
use a combination of Mean Field (MF) analysis and Monte Carlo simulations to characterize the 
density profile in growing lattice and construct a MF phase diagram. This study will have possible 
implications for transport and patterning in growing fungal filament. 



An intrinsic feature of the biological world is the process of growth. One fascinating aspect of growth process 
during embryonic development is the formation of patterns. Positional information associated with pattern formation 
is initiated by concentration gradient of morphogens [TJ |2] . Experimental studies such as the one done with wild-type 
D. melanog aster embryos have indicated that the bicoid morphogen distribution scales linearly with the size of the 
embryo [2 , which would imply that as a function of the relative position in the embryo, the morphogen distribution 
is invariant. In general such scaled patterning is seen to be robust against environmental variations and genetic 
background variations [2j. A natural question that arises in this context is - what are the mechanisms by which 
such scale invariant patterning is achieved ? The answer to this question is an involved one, which requires careful 
consideration depending on the specific biological context. While we are not going to address that issue in this letter, 
it would serve as the broad motivation for focusing our attention on a simple one dimensional lattice model for growth 
which exhibits such scaled growth pattern. For such scaled growth pattern, the density profile remains invariant in 
the scaled position variable x/L, where x is the position and L the length of the lattice. This model that we would 
discuss, incorporates some of the ingredients involved in quasi one-dimensional biological growth process of individual 
fungi filament. 

Fungi are eukaryotes and can typically grow as cylindrical tubes (hyphae), with the growth happening at the tip of 
fungi [3]. The material required for growth is packaged as vesicles that are supplied to the growing tip by molecular 
motors. The molecular motors transport these materials using the microtubules tracks to the growing end of the 
hyphae. These vesicles bring fresh membrane, and regulators neccesary for growth. The growth process of fungal 
hyphae have a sense of persistence in one direction. In addition to linear growth of the filament, lateral branching of 
hypha also appears giving rise to a network termed as mycelium [3]. Fungal mycelium exhibits patterning and it has 
also been observed that the spatial pattern associated with morphogenesis in fungi during spore germination requires 
the presence of microtubules [4 . 

We model the growing fungal filament as a one dimensional lattice which can grow at the one of the ends of the 
lattice. The model that we study belongs to the class of driven lattice gas models and a generalization of the model of 
dynamic lattice first proposed and studied in [5j [6] . Driven lattice gas models have been used to model a number of 
biological processes, such as the motion of ribosomes in m-RNA |7 , motor and vesicle transport [5HT2]. transport across 
biomembranes [13 and growth of fungal hyphae [SJ|6]. Many of these models are based on the totally asymmetric 
exclusion process (TASEP), in which the particles jumps unidirectionally on the lattice with a single rate and interacts 
with other particle via hard-core repulsion. For TASEP, the entire non-equilibrium phase diagram in terms of particle 
fluxes at the boundaries, has been worked out exactly [14]. TASEP has been generalized to incorporate a dynamically 
extending lattice, which is growing at one the ends of the lattice [SJIS]- For such dynamically extending TASEP, the 
growth velocity at tip is proportional to the particle density at the tip and this minimal model has been proposed 
to describe the growth of fungal hyphae [5j [6] . Number of other kinds of dynamic lattice gas models have studied in 
varied context such as a dynamic lattice coupled to a diffusive reservoir [T5HT9] and a dynamic lattice which includes 
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process of corrosion [20]. 

Here we consider the coupling of the lattice growth process with Langmuir kinetics. For the case of growth of 
fungal hyphae, the fact that the transported cargo may detach (and reattach) to the microtubule filament makes it 
plausible scenario to explore. In particular, we study the dynamic regime where these two processes compete; the 
kinetic rates are such that fluxes of the particles at the boundaries are comparable to that fluxes in the bulk. For a 
static lattice, the competition beween particle fluxes at the boundaries with the process of particle exchange with the 
surrounding bath has been addressed in [2T] , We study this competition in context of dynamically extending lattice, 
focussing on a regime where the site attachment and detachment rates scale inversely as the system size. We use mean 
field (MF) analysis to study the system and compare these results with those obtained by perfoming Monte Carlo 
simulations. We find that the competition between the growth process and kinetics of attachment and detachment 
not only leads to the possibility of phase coexistence and shocks in the co- moving frame of the growing lattice, it also 
implies scale invariant density profile in the growing lattice. This dynamic regime exhibits qualitatively new scenario 
when compared to the situation when the growing lattice is coupled to a bulk reservoir where the total attachment 
and detachment rates scales extensively with the filament size [T7] . 
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FIG. 1: Schematic representation of the dynamical processes on the growing lattice with particle entry rate a, growth rate 7, 
attachment rate ujd, detachment rate uo e and hopping rate 1. 



We represent the growing filament at any instant of time to be a one-dimensional lattice of length L and lattice 
sites N with a lattice spacing e = L/N. Particles hop unidirectionally to the left on this lattice with the instantaneous 
configuration being described in terms of the occupation number of particles at each lattice site z, which can assume 
values rti = 1 if the site i is occupied or = if the site is vacant. The lattice length and the number of lattice sites 
increase due the growth process which happens at the left end of the filament. The lattice growth happens when a 
particle moving unidirectionally towards the left on the lattice, reaches the filament end and creates a new lattice site 
with a rate 7 so that the total lattice length increases by an amount e. The non- growing right end of the filament is 
supplied with particles which enters the rightmost site with a rate a. In the bulk of the lattice two distinct processes 
can occur: particles can hop to the neighboring site on the left with rate 1, if the site is vacant; particles can detach 
from the filament with rate uo e or a new particle can attach to the filament lattice site with a rate 00^ provided the 
site is vacant. The various dynamical processes are summarized in Fig.l. We focus on the situation, similar to [21], 
where the overall (un)binding rates for the entire lattice is held fixed, so that the site (un)binding rates vary inversely 
with the system size. Since the lattice size continously increases, there are two frames to consider; while one is the 
laboratory frame, the other is a frame attached to the growing tip. We choose the reference frame attached to the 
growing tip of the filament, with the leftmost site being labelled as 1 and the rightmost being labelled as site N. We 
write the evolution equation of the occupation number of particles in the aforementioned frame. Any site i measures 
the distance of that site from the growing tip. Each lattice site can be occupied by at most one particle. Whenever 
the lattice grows due to creation of a new site at the tip, all the site labels are updated i — >> i + 1. The evolution 
equation for the occupation numbers in the bulk at site i is, 

= 7l i+ i(l -7li) -7li(l -Tli-x) + 777,1 (7li_i -m) 

- u e ni + u d (l-rii), i>3 (1) 

The terms on the right hand side of the equation have their usual interpretation of gain and loss terms due to 
translation of the particles, the effect of the growing tip and the Langmuir kinetics. The corresponding current at 
site i can be identified as, Ji = (1 — rii-\)ni — n\Ui-\. The occupation number of the particles at sites 1 and 2 in the 
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comoving frame of the tip are respectively. 



dt -7ni +n 2 (l -ni) (2) 
= n 3 (l - n 2 ) - n 2 (l - ni) - <ynin 2 - u e n 2 + u d (l - n 2 ) 



At the right end of the filament the evolution equation reads, 



a(l - n N ) - n N (l - un-i) (3) 

Defining p = (n^), we introduce mean field approximation which amounts to taking averages over ensembles and 
factorizing the two-point correlators, (n^+i) = (n^)(n^ + i) [121 [22]. We now write the continuum evolution equation 
in terms of the rescaled variable x = ie/L, which measures the relative position of the particle in the growing lattice. 
We define the reduced rates, Qd = Nuod and Vt e = Nuo e . These rates have the interpretation of the total bulk 
attachment and detachment rates of particles for the entire growing lattice. In the thermodynamic limit of N — >• oo, 
we obtain the the continuum evolution equation as, 



%=iw-p)\-^%-^ P +w-p) (4) 

where we have retained terms upto leading order in 1/N. 

In the continuum description, the evolution equation of density at the non-growing end at x = 1 satisfies, 

^1 = a[l-p(l)]-p(l)[l-p(l)] (5) 
The evolution equation for particle densities at sites 1 and 2 are respectively, 



dp 

~dt 
dn 2 



Of = -7Pi +P2(l-Pi) (6) 



df = P3(l - P2j -P2(l -Pi) -7P1P2 (7) 

where p\ = (m) and p 2 = (n 2 ). In order to determine the steady state density at the boundary of the growing tip, 
we make the approximation, 



Then using Eq.([6| and Eq.0 we obtain, 



Ps = P2 (8) 



P2 = I-27 (10) 
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The approximation of Eq.([8| works reasonably well, when the bulk phase adjoining the tip of the growing tip is a 
high density phase. So for example if the bulk phase is a high density phase or even a mixed phase of high and low 
density, then this approximation can predict the steady state density profile at the growing tip boundary reasonably 
well, as seen by comparison with profiles obtained by Monte Carlo simulations. We shall elaborate more on the 
domain of validity of this approximation, when we discuss the possible phases and construct the phase diagram later. 

Having obtained the the expression for pi, the continuum steady state equation assumes the form, 



dJ 

dx 



where, 



(n d + n e ) P + n d = 



1-7 

J has the interpretation of particle current in the bulk of the lattice in the frame of the moving tip. 



(ii) 



(12) 



Integrating Eq.(ll), we obtain, 

2(P6 - p) + 



7(1 - 27) 2n d 
1-7 n 



log 



Vtp h - Q d 



Q(x — Xb) 



(13) 



where p^ is the boundary density and £5 is the position of the boundary; £5 = for the left boundary and £5 = 1 
for the boundary at the right end of the filament. Here Q — Q d -\- Q e . 
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FIG. 2: Steady state density profile of particles in scaled variable x. Here, a = 0.15,7 = 0.1, fid = 0-2 and fi e = 0.3. The 
square, triangle and circles with solid lines represent the density profile obtained by Monte Carlo numerical simulations for 
different sizes of the growing lattice. The dashed lines are the MF densities solutions corresponding to low density(LD) and 
high density (HD) phases obtained from Eq. (13). 



When the density profile satisfies the boundary condition at x = 0, then the boundary density is obtained by 
equating to the occupation number density at site 2 to the boundary density, so that p(0) = p2 = 1 — 27. This density 
profile corresponds to a High density (HD) phase. For the density profile which satisfies the boundary condition at the 
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other end of the filament located at x = 1, we obtain p(l) = a. This density profile corresponds to a Low density (LD) 
phase. In general, the boundary conditions imposed at the two filament ends cannot be satisfied simultaneously. As 
a result, at least two domains usually develop along the filament, giving rise to the possibility of shock profiles. The 
location of the shock is determined by equating the current J for the LD and HD solution at the position of the shock. 
Fig. [2] shows shock density profile obtained by Monte Carlo simulations for different system sizes. These profiles are 
compared with the MF prediction of Eq.(13) corresponding to the LD and HD mean field solution of density. As seen 
in Figj2j the entire density profile including the location of shocks overlaps for the different system size in the scaled 
variable x. This would imply that the growth process is such that the density profile remains invariant for different 
system sizes in the scaled variable x. We also observe that the MF predictions matches reasonably well with the 
Monte Carlo simulation results. In order to analyse the different dynamical growing phases that are possible for this 
system, we focus on a special regime where the kinetic rates satisfy the constraint, 



2 £l e + £l d 



K d 



For this choice of constraint, the steady state equation for the bulk reduces to, 



d 



dp 

dx 



B 



= 



(14) 



(15) 



where B = (Q d + Q e )/2. 

Thus for this particular choice of constraint, the steady state MF density profiles either vary linearly with x or are 
homogenous with respect to x. As a result, there are three different profiles which may coexist along the filament. The 
different solutions correspond to a linearly decreasing high density (HD) solution which satisfies the boundary density 
at the growing tip (x = 0), a linearly decreasing low density (LD) solution which satisfies the boundary condition at 
x — 1 and a homogenous density solution (U) which is similar to a density profile determined by Langmuir kinetics 
alone. The non-equilibrium phase diagram can be constructed once we determine the limits of coexistence of the 
phases. Phase coexistence in bulk can occur when the densities in the linear phases (LD or HD) matches with the 
density of the homogenous solution (U) at x — x q , m the domain range < x < 1, where x Q is position of the domain 
wall separating the two phases. For LD phase to coexist with the HD phase, current of the LD phase must equal the 
particle current in the HD phase at domain wall position x — x Q in the domain < x < 1. Phases on the either side 
of the domain wall would have a density profile corresponding to the particular phase on the growing lattice. The 
phase boundaries are determined by the condition of expulsion of one of the phases from the bulk, i.e; by setting the 
location of the domain wall position to the filament boundaries. While the density changes continuously on moving 
across the phase boundary separating the U region with the HD(LD) region, the density change is dicontinuous while 
moving from a LD region to a HD region. 



In order to analytically determine the MF density profile and the phase boundaries, separating the different phase 
regions, we again use the approximation introduced in Eq.([8|, so that the expressions for pi and p2 are determined by 
Eq.([9| and Eq.(10), as before. This approximation scheme implicitly assumes that pi is determined by the boundary 
conditions at the growing end of the filament alone. This is valid when the bulk density adjoining x = is determined 
by the growth rate at growing top alone, and the rate of tip growth is the rate limiting process which controls the bulk 
density adjoining the region of the growing tip. However when entire bulk phase is a LD phase and the HD phase 
is totally expelled from the growing end of the filament, then this approximation will not hold good. As long as the 
high density phase is present in the bulk, this approximation is reasonable, as confirmed by Monte Carlo simulations 
and illustrated in Fig. 2. The condition for linear solution is then, 



7(l-2 7 ) " 
1-7 _ 



(16) 



and the corresponding expressions for MF density profile for the LD, U and HD phases are respectively, 
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p[x) = < 



' Phd(x) = 1 -2~f + Bx 
o (x) - 1 - 7(1 ~ 27) 

AM X J — 2 2(1-7) 

, Pid(x) = a + B(l-x) 



(17) 



By varying the parameters 7 and a, for a fixed value of f^, one can see the existence of both single phase regions, 
regions of phase coexistence and presence of shocks in the bulk. 



Phase coexistence line between the HD and LD-HD: This is determined by using the expression for current J in Eq.([12| 
and matching the current for the HD phase, with the current in LD phase at x = 1, so that J^(l) = 7/^(1). We thus 
obtain the equation of phase boundary as, 



where, J hd (l) = (1 - 2 7 - B)(2 7 + B) 



2 I-27 + 27- 
1-7 

7 (l-2 7 )(l-2 7 -B) 
1-7 



OL + JffD(l) =0 



(18) 



Phase coexistence line between the LD and LD-HD: On matching the current for the HD phase, with the current in 
LD phase at x = 0, so that J^(0) = J/i</(0), we obtain the equation of phase boundary as, 



a 



1-2B 



7(1 - 2 7 ) 



1-7 



a + J hd (0) + B 2 -B 



1 - 27 + 2 7 2 
1^7 







(19) 



where, J^(0) 



7(1-27) 
1-7 • 



Here we have determined the location of phase boundaries in the phase diagram using 
Eq.([16j). This maybe understood in the following manner: Consider a region in the vicinity of the left boundary with 



x = S. Let the region on the right of x be a LD phase while the region on the left (from x = to x = S) be in the HD 
phase. For this situation, condition for the phase boundary separating the LD phase with LD-HD phase coexistence 
region is obtained by letting the domain wall position x Q — >• 0. For a finite value of 5, p\ = \^ and the condition of 
linearity for the coexisting LD and HD profile is set by it. Consequently if one approaches S — >• from the right, the 
boundary density for site 1 and 2 can be obtained using Eq.([9| and Eq.(10) and the condition for linear profile is set 



by Eq.(16). Since the phase boundary can be obtained by the limiting procedure of setting the domain wall position 
x Q to x = 0, approaching from the right, thus the phase boundary itself may be determined using the same condition 
of linearity as Eq.(16) and value of p\ is determined by eq.d9|) although the LD phase on the otherside of the phase 



boundary in the phase diagram would not satisfy the same condition for linear solution as before. 



Phase coexistence line between the LD-HD and LD-U-HD: If we require that the density for LD phase is same as 
that of U phase at the domain wall location, then the expression for position of the domain wall is xid/ u = a ~^ B B ~ pu . 
Similarly expression for the domain wall position separating HD with U region is Xhd/u — 1 ~ 2l B ~ pu . The condition 
for phase coexistence boundary is that x\^i u = x^d/u- Thus we obtain the equation for phase boundary as, 



a = 1 - 27 



^d 
2 



"2 - 2 7 -7 2 " 
_1 - 2 7 + 7 2 _ 



(20) 



Phase coexistence line between the HD and HD-U: By matching the density for the HD phase with the density in U 
phase at x = 1, we obtain the equation in terms of 7, 



2 7 2 + (4 - 2^)7 + 1 - 2B = 



(21) 



7 



Since B = Qd + Qe = ^2^+272 ■> we numerically solve Eq.(|2l|) to determine the value of 7 corresponding to the phase 
boundary. 

Phase coexistence line between the U and HD-U: This is determined by matching the density for the HD phase, with 
the density in U phase at x = 0. Thus the equation of phase boundary is, 



Thus we see that above this critical value of 7 determined by Eq.(22), the high density phase is expelled from the 
bulk of the lattice. In such situation the only possibilities are that of homogenous U phase, a pure LD phase and 
LD — U phase coexistence. 

For a pure LD phase, the rate limiting process is the input rate a, so that the linear profile, p(x) = a — B(l — x) 
would occupy the entire bulk of the lattice, right upto the growing end of the filament. In such circumstances, in 
order to determine the boundary value of density at x = 0, we use the approximation, 



p 2 =a+^(n d + n e ) = a (23) 

This approximation also would imply that p\ = . Thus the corresponding linear solution for the low density 
profile would then have to satisfy the condition, 



K d =Ul--^-) (24) 
2 V j + a Q J 



This condition for linear solution is different from the condition of linearity expressed in Eq.(16), where the HD phase 
was present in the bulk. 

Phase coexistence line between the LD and LD-U: This is determined by matching the density for the LD phase, 
with the density in U phase at x = 0. The value of density for the LD mean field solution at x = is, pu(0) = 
a + -^(fLd + O e ). The homogenous solution which determines the U phase density is given by p u (0) = |(1 — 7Pi). We 
use the approximation of Eq.(23) and the steady state condition corresponding to Eq.(|6]), to obtain the corresponding 



expression for the p u (0). By equating the expressions for the densities in the two phases, we obtain we density at the 



left boundary as p = |[1 — 37 + \/~f 2 + 27 + 1], and the equation for the phase coexistence boundary is, 

a = p (25) 

7 + 9 ~ 19 

The corresponding condition for linear solution is, a = — 2 n^K d • 

Phase coexistence line between the U and LD- U: The coexistence boundary is determined by matching the density for 
the LD phase with the density in U phase at x = 1. For this case 9id(X) = ct and p u (l) = \(1 — jpi) = \{1 — ^q^), 
where we have again made use of the approximation, 

P2 = a (26) 



Thus the equation of phase boundary is, 

2a 2 — a 
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FIG. 3: Phase plane cut along the a — 7 plane on varying the overall attachment rate of particles. (a)Qd = 0: This corresponds 
to the situation where Langmuir kinetics is switched off. The dashed curves are the MF results derived in Ref.[6 . (b) Qd = 0.1: 
Regions of phase coexistence are seen, (c) Qd = 0.3 : The LD-U region widens towards the left while HD-U region expands 
downwards and LD phase region is almost expelled, (d) Qd = 0.5 : HD phase is expelled from the phase diagram. The location 
of (U)-(LD-U) and (U)-(HD-U) phase coexistence line remains unchanged on variation of Qd in all cases. 



The corresponding condition for linear solution is, a = The MF phase diagram in Fig. [3] displays the phase space 
cut along a — 7 plane for a special choice of parameter values which are consistent with the constraint requirement 
specified by Eq.( 14). Fig j3] illustrates how with the increase of regions of phase coexistence in the bulk is expanded. 
When the Langmuir kinetic process is switched off, then the dynamics of the model reduces to the case of a dynamic 
TASEP studied in Ref.©. Fig. ga) compares the phase boundaries obtained for the case when Vt^^e — with 
the results obtained in Ref.[6 for dynamic TASEP . While the HD-U and LD-HD phase boundary obtained in this 
case is exactly the same as that obtained in Ref.[6 , the location of LD-U phase coexistence lines differ due to the 
approximations made through Eq.(26). Although the aspects of phase coexistence and presence of density shocks in 
the bulk, have many similarities with the phase diagram for TASEP-LK process for static lattice [8 , the topology 
of the phase diagram in this case differs significantly. Comparison of the results obtained by this approximate MF 
analysis with Monte Carlo simulations reveal that while their is broad qualitative agreement for the density profiles 
and the rough location of the phase boundaries, the MF density profile does not exactly match with the actual 
density profile obtained by Monte Carlo simulation. The MF analysis is also not able to accurately determine the 
exact location of the phase boundaries. However the feature of scale invariant density profile predicted by the MF 
analysis is robust as seen by comparison with Monte Carlo simulations done for different system sizes. 



In summary, we have presented a model of dynamically extended lattice which incorporates the process of Langmuir 
kinetics and exhibits feature of scale invariant density profile in the bulk. It also rationalizes in simple terms the 
interplay of the various competing processes which affect growth. In particular it highlights the dynamic regime for 
which the bulk process of Langmuir kinetics competes with the boundary processes of particle input and growth at 
the tip. The resultant phase diagram shows rich behavior, which includes different inhomogeneous phases, phase 
coexistence in bulk and the possibility of density shocks. This model will have significant ramification in modeling 
transport and growth in fungal hyphae. The emergence of scale invariant patterned growth which can be altered by 
the interplay of internal regulatory mechanisms with environment effects, makes it a useful minimal model which 
may be suitably modified for modeling branching and patterning in fungal mycelium. The framework that we discuss 
provides a natural means to link internal regulatory mechanism like vesicle supply rate with the branching or tip growth 
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rates which can be dependent on external cues [3 j. The possible extensions of this model can involve extending the 
model to include bidirectional transport along the filament and a more realistic biological scenario would also involve 
generalizing it to account for transport on multiple parallel tracks, as is the case for transport in fungal hyphae [3]. 
It would then be interesting to study whether transport on multiple tracks can lead to phase segregation of material 
within the tracks as has been seen in [23] . 



I. ACKNOWLEDGMENTS 



I would like to thank G. Tripathy for useful discussions and suggestions. 



[1] N. H. Patel and S. Lall; Nature 415, 748 ( 2002) 

[2] D. Cheung, C. Miles, M. Kreitman and J. Ma, Development 138, 2741 ( 2011) 

[3] J. Deacon, Fungal Biology (Blackwell, Oxford, 2006, 4th ed) 

[4] S. D. Harris, Mycologia 97, 880 ( 2005) 

[5] K. E. P. Sugden, M. R. Evans, W. C. K. Poon and N. D. Read, Phys. Rev. E 75, 031909 (2007) 

[6] K. E. P. Sugden and M. R. Evans, J. Stat. Mech. P11013 (2007) 

[7] C. T. MacDonald, J. H. Gibbs J and A. C. Pipkin, Biopolymers 6 1 (1968) 

[8] A. Parmeggiani, T. Franosch and E. Frey, Phys. Rev. E. 70 046101 (2004) 

[9] Y. Aghababaie, G. I. Menon and M. Plischke, Phys. Rev. E. 59 2578 (1999) 

[10] P. Ashwin, C. Lin and G. Steinberg, Phys. Rev. E. 82 051907 (2010) 

[11] S Muhuri and I. Pagonabarraga , EPL. 84 58009 (2008) 

[12] S. Muhuri and I. Pagonabarraga, Phys. Rev. E. 82 021925 (2010) 

[13] T. Chou and D. Lohse, Phys. Rev. Lett. 82 3552 (1999) 

[14] M. R. Evans, D. P. Foster, C. Godreche and D. Mukamel, Phys. Rev. Lett. 74 208 (1999) 

[15] M. Ebbinghaus and L. Santen, J. Stat. Mech. P03030 (2009) 

[16] M. Ebbinghaus, C. Appert-Rolland and L. Santen, Phys. Rev. E. 82 040901 (2010) 

[17] K. E. P. Sudgen, Ph. D Thesis, (University of Edinburg, U.K, 2009) 

[18] S. A. Nowak , P. W. Fok and T. Chou, Phys. Rev. E. 76, 031135 (2007) 

[19] J. Tailleur, M. R. Evans and Y. Kafri, Phys. Rev. Lett. 102, 118109 (2009) 

[20] S. Dorosz, S. Mukherjee and T. Platini, Phys. Rev. E 81, 042101 (2010) 

[21] Parmegianni A, Franosch T and Frey E, 2003 Phys. Rev. Lett. 90 086601 

[22] Non equilibrium Statistical Mechanics in One Dimension, edited by V. Privman ( Cambridge University Press, Cambridge, 
1997) 

[23] S. Muhuri and I. Pagonabarraga, J. Stat. Mech. P11011 ( 2011) 



